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Effects of possible orbital order in magnetic properties of two-dimensional spin gap system for 
CaV409 are investigated theoretically. After analyzing experimental data, we show that single 
orbital models assumed in the literature are insufhcient to reproduce the data. To understand 
j/j ' the origin of the discrepancy, we assume that in d^ state of V, dxz and dyz orbitals have 

substantial contributions in the lowest-energy atomic level which leads to a double-degeneracy. 
We study possible configurations of the orbital order. By exact diagonalization and perturbation 
expansion, we calculate the susceptibility, wavenumber dependence of low-lying excitations and 
^ , equal-time spin-spin correlations which is related to integrated intensity of the neutron inelastic 

O \ scattering. These quantities sensitively depend on the configuration of the orbital order. The 

calculated results for some configurations of the orbital order reproduce many experimental 
results much better than the previous single-orbital models. However some discrepancy still 
^ ' remains to completely reproduce all of the reported experimental results. To understand the 

.^ ■ origin of these discrepancies, we point out the possible importance of the partially occupied dxy 

^D ■ orbital in addition to orbital order of partially filled dxz and dyz orbitals. 
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'^ ' §1. Introduction 

O ■ Since the spin gap was discovered in a typical two-dimensional material CaV409,El^ the mechanism 

of the spin gap formation in this system has been intensively studied by many theoretical methods. 
ii@ii0ii0§) 



This material consists of VO5 pyramid layers. Crystal structure of the VO5 pyramid layer is 
shown in Fig. 1(a). The oxygens constitute a complete square lattice while in the square lattice 
of the vanadium atoms, 1/5 of them are depleted. The vanadium atom is nearly located at the 
center of the pyramid constructed from four oxygens in the layer and one apical oxygen. Each 
VO5 pyramid is connected by edge sharing. The unit cell of the layer includes two edge-shared 
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plaquettes of V atoms shown in Fig. 1(a), because the apical oxygen of one edge-shared plaquette 
is located above the VO2 plane while that of the other one is below the plane. 



Fig. 1. (a) Crystal structure of VO5 pyramid layer. Full circles represent vanadium atoms, while open circles are 
oxygen atoms. Apical oxygens are omitted in this figure. Shaded square area represents the unit cell of this 
system. Vectors ai and 61 are the unit lattice vectors ; ai — 3x — y and 61 = a; + 3y, respectively. Here x and 
y represent the unit lattice vectors of the square lattice composed of oxygens defined as x = ae^ and y — aSy. 
Parameter a is a lattice constant while e^ and By are unit vectors, (b) Four kinds of the superexchange couplings. 
Bold-gray, bold-black, thin-gray and thin-black lines represent the spin exchange couplings, Jep, Jg^j, Jcp and J^j , 
respectively. 



Since the valence of V atom is 4-|-, the d electron on V atom can be treated as a nearly localized 
spin with S = 1/2. In the literature, it has been assumed that a non-degenerate orbital on V atom is 
occupied.oB'Q'U'Q'Ll'U'EI'llj'lllP Based on this assumption, the S = 1/2 antiferromagnetic Heisenberg 
(AFH) model with the nearest-neighbor and the next-nearest-neighbor exchange couplings has 
been introduced. When the nearest-neighbor exchange couplings are only taken into account, 
the system is described by the AFH model on the square lattice of the plaquettes, where each 
plaquette consists of four V atoms almost centered in the edge-shared pyramid. (We call this lattice 
plaquette lattice.) On the other hand, when only the next-nearest-neighbor exchange couplings are 
considered, the system is described by two disconnected square lattice of the plaquettes, where 
each plaquette is twice as large as and 45° tilted from the plaquette element of the previous case. 
In this case each plaquette consists of four corner-shared pyramid. Here, four different kinds of 
spin exchange couplings are introduced, as is shown in Fig. 1(b). The spin exchange couplings 
between V atoms in the edge-shared and the corner-shared plaquettes are represented by Jep and 
Jcp, respectively. While the spin exchange couplings connecting the edge-shared and the corner- 
shared plaquettes are J^^ and J^^, respectively. In the literature, the case where Jg^^ ~ Jep 
is taken as J and J^^j ~ Jcp is J', has been mainly studied theoretically since only two kinds 
of spin exchange couplings should be derived from the previous single-orbital assumptions. The 
effects of the Jahn- Teller distortion, tilting of the pyramids and the spin-orbit couplings may 
cause the differences between Jg^^ and Jep or between Jg^^ and Jcp. The theoretical methods such 
as perturbation expansionspHo' exact diagonalization,EP mean field approximations,LrcP quantum 
Monte Carlo methodpEP high-temperature expansionsu'Q'B^ and the density matrix renormalization 
group methodE3 have led to the result that this system has a spin gap in the region around 
Jep > Je(j without frustration (J' = 0) or in the region J' ~ 0.5J. Here the spin gap is defined as 
the energy difference between the singlet ground state and the lowest triplet states. In the above 



mentioned regions, the origin of the spin gap is ascribed to the edge-shared plaquette singlet. The 
energy dispersion of the triplet states has been calculated by the perturbation expansion from the 
edge-shared plaquette singletfloQ' The wavenumber of the lowest energy excitation is (vr,7r) in 
the absence of the frustration, which is also obtained by the variational Monte Carlo method.O 
However, it shifts to incommensurate wavenumbers as the frustration due to J' becomes relatively 
large. 

Recently, the neutron inelastic scattering study for single crystal of CaV409 has been performed. 
113 i3 The experiments indicate that the predictions from the theoretical studies for the conventional 
model are inconsistent with the experimental results. First of all, the wavenumber of the lowest 
energy excitation given by the experiments is (0,0) in the magnetic first Brillouin zone, which is 
different from the theoretical results. Here, the magnetic first Brillouin zone in the wavenumber 
space is expanded by the vectors (12 = di + bi and 62 = —di + ^ii as shown in Fig. 2. It is \/2 x \/2 
times larger than the first Brillouin zone of the unit cell expanded by the reciprocal lattice vectors 
tti = ai/lOa^ and 61 = bi/lOa^, where a is a distance between the nearest-neighbor V atoms. 
Secondly, if the spin exchange couplings between V atoms are determined so as to reproduce the 
dispersion of the lowest triplet excitations in the experimental results, a large difference is resulted 
between the couplings, i.e. J^.^^ = 0.088Jcp and Jep = Jq^ = 0.395Jcp,c3) which is hard to 
understand within the conventional framework. 



Fig. 2. Brillouin zone of two-dimensional VO2 plane in CaV409 lattice. Square surrounded by broken lines is 
the first Brillouin zone of the lattice, while the square with bold-solid lines is the magnetic first Brillouin zone 
obtained by the experiment .ta' The wavenumber q is represented by the vectors hi and bi. While fc is described 
by £i2 = di + bi and 62 = — di + bi, respectively. Here, the lattice constant between the nearest-neighbor V atoms 



In order to understand the contradiction between the experimental results and the theoretical 
ones, we should reconsider the mechanism of the spin gap formation for CaV409. In this paper, we 
study effects of orbital degeneracy and orbital order which has been neglected in the literature. The 
crystal field from the oxygen ions on the corners of the pyramid lifts the degeneracy of t2g orbitals in 
the atomic level. However, two orbitals whose wavefunctions are expanded by d^z and dy^ orbitals 
may still be degenerate even in this crystal field. Then the effective Hamiltonian sensitively depends 
on configuration of the occupied d orbitals through the mechanism of the superexchange interaction 
if the occupied orbitals have substantial contributions from the dxz and dyz orbitals. In Section 2, 
we introduce several possible effective spin Hamiltonians with the orbital order. It is found that 
strength of spin exchange couplings strongly depends on patterns of the orbital occupancy. This 



may explain the appearance of large difference in the spin exchange couplings needed to reproduce 
experimental results. In Section 3, we estimate the strength of the spin exchange couplings from 
analyses of the temperature dependence of the uniform magnetic susceptibility. We use the exact 
diagonalization (ED) method. Although the system size tractable by the ED is limited, the cases 
we studied show rather small system size dependence due to the opening the spin gap which makes 
it possible to infer the thermodynamic limit. Then we calculate the energy dispersion of the triplet 
states and the wavenumber dependence of the equal-time spin-spin correlation which is related to 
the integrated intensity of the neutron inelastic scattering within the perturbation expansion (PE). 
The important elements to explain the experimental results are pointed out. The origin of the spin 
gap in each case of the orbital order is also discussed. In Section 4, we compare these theoretical 
results with the experimental ones and discuss the importance of the orbital order. Section 5 is 
devoted to summary. 

§2. Eflfective Hamiltonians 

For the purpose of understanding the problems mentioned in the previous section, we consider 
several models for the spin gap formation. The Hamiltonian is given by the AFH model written as 

H= J2 JijSvSj, (2.1) 

where Si represents the spin operator with S = 1/2 at i site and Jij is the spin exchange coupling 
between the spins at i and j sites. 

In this study, we assume that the orbitals giving the lowest energy in the atomic level due to the 
crystal field are doubly degenerate. In other words, the ground state wavefunction has substantial 
weight of the dxz and dyz orbitals and that the d electron occupies these orbitals at least partially. 
Below we refer to occupied orbitals simply as dxz or dyz orbitals in this paper. However, it does 
not necessarily mean that those consist of pure dxz and dyz orbitals. What we need to make the 
proposals in this paper relevant is the double-degeneracy of occupied orbitals with dxz and dyz 
component contained. Even small contribution of dxz and dyz orbitals may make the exchange 
coupling sensitively dependent on the occupied orbitals. Although the static Jahn- Teller distortion 
or tilting of the pyramid could lift the degeneracy of dxz and dyz orbitals, we do not consider this 
effect within this study, because such static distortions have not been observed so far. Therefore, 
we investigate effects of the orbital order for the dxz and dyz orbitals. As is mentioned later, the 
strength of the spin exchange couplings between V atoms depends on the configuration of the 
occupied d orbitals on V atoms. 

We consider the spin exchange couplings between the nearest as well as between the next-nearest 
neighbor V atoms through the superexchange mechanism, as is shown in Fig. 3. In the first case, the 
V-O-V bonds make approximately right angle. Then the superexchange coupling works relevantly 
through the Pz orbital of the oxygen when the d electron occupies dxz orbital on V atom and 



Fig. 3. Lattice structure of V-O-V and configuration of occupied orbital. Figures (a) and (b) show the cases with 
the right angle of V-O-V. Figure (a) shows the case that the occupied d orbitals on V atoms are different, while (b) 
is the case with the same occupied orbital. Figures (c) and (d) show the cases with the straight line of V-O-V. The 
configurations of the occupied d orbitals in Fig.(c) and (d) are the same as those in Fig. (a) and (b), respectively. 



the other d electron is locaHzed at dyz orbital on the other V atom shown in Fig. 3 (a). However, 
the superexchange coupling may be relatively small when the both d electrons occupy the same 
orbitals on each V atom shown in Fig. 3(b). The reason is that one of the transfer integrals between 
Pz orbital on oxygen and one of the d orbitals is relatively small due to the symmetry of the 
wavefunction of these orbitals. In the second case, the V-O-V bonds make approximately straight 
line. Then the superexchange coupling works effectively when the d electrons are localized at the 
same orbitals which extend to the oxygen as shown in Fig. 3(d). While the superexchange coupling 
should be relatively small in the case shown in Fig. 3(c). When the d electrons occupy the orbitals 
both of which do not extend to the oxygen located between the V atoms, the strength of the 
superexchange couplings should be also relatively small. Consequently, the anisotropy of the spin 
exchange couplings arises through this mechanism. 

Here, we introduce five models with the possible patterns of the orbital order as shown in Figs. 4. 
The orbital patterns in Figs. 4(a) and (b) are the cases where the electrons occupy the d orbitals 
alternatingly within the edge-shared plaquettes but in a uniform way from plaquette to plaquette. 
The relevant spin exchange couplings in Fig. 4(a) are the edge-shared plaquette bonds Jep, while the 
ones in Fig. 4(b) are Jcp and J^^. Figure 4(c) shows the case where the plaquette unit in Fig. 4(a) 
and (b) are placed alternatingly. Figure 4(d) shows the case where all the occupied orbitals have 
the same symmetry. Figure 4(e) shows the case where the occupied d orbitals are alternating in the 
x-direction while are uniform in the ^/-direction. From here, we take the strengths of the relatively 
small spin exchange couplings as the same values. In reality, these value are not necessarily the 
same because these small couplings appear due to some small distortion of the lattice from the 
perfect square lattice. Then to study difference between them, these small effects have to be 
considered. This point is discussed in detail in Section 4. In other possible patterns of the orbital 
order, the magnetic unit cell becomes more complicated. In this paper, we do not consider these 
more complicated possibilities. 

§3. Results 

In this section, we study magnetic properties for the models in Figs. 4(a), (b) and (c) in more 
detail than for (d) and (e) since the configuration of the occupied orbitals in the edge-shared 



Fig. 4. Five possible patterns of orbital order viewed in the projection to the xy plane. Bold lines connecting V 
atoms represent the dominant superexchange couplings. The symbols on V atoms represent the dxz and dy^ orbitals 
schematically, (a) A case with the alternating pattern of the orbitals in the edge-shared plaquette. The filled orbitals 
are "tiled" in a uniformed way among different plaquettes. (b) This is the same as the one in Fig. (a). However, the 
dxz and dyz orbitals are replaced each other, (c) The case where the orbital in all the nearest-neighbor V atoms 
are occupied always in the alternating configuration, (d) The case that all the orbitals are occupied uniformly, (e) 
The case that the occupied orbitals are alternating in the a;-direction while are uniform in the y-direction. 



plaquette for these models has four-fold rotational symmetry which appears to be consistent with 
experimental observations .I12P 

3.1 Uniform Magnetic Susceptibility 

In this subsection, we determine the amplitude of the relevant spin exchange couplings and 
the other relatively small ones to reproduce the temperature dependence of the uniform magnetic 
susceptibility x given by the experiments. Here, we assume that eq.(p.l|) is an effective spin Hamil- 
tonian of magnetic properties below 700K. We calculate the x of 16-site system by the ED with 
the periodic boundary condition. The uniform magnetic susceptibility is defined as 

^ /3 TVE.E,-^fg|exp(-/?6„) 
^ N Trexp(-/3e„) ' ^ ' ^ 

where f3 represents the inverse temperature and A^ is the number of site. In terms of comparison with 
the experimental data in the unit of emu/g, we should multiply our data with energy scale J by a 
factor A:Np^g^^-Q/ Jk-QM with N^ the Avogadro number, /ig the Bohr magneton, fcg the Boltzmann 
constant and M the gram per mole. The factor 4 comes from the number of V atoms in a unit cell. 
The parameter J depends on each model. Here, the g-ia.ctoi is taken as a fitting parameter. The 
experimental data show a peak at a temperature about Tp ~ llOK. The temperature T* where the 
amplitude of x becomes a half of that at Tp is about 595K. We choose the amplitude of the spin 
exchange couplings so as to give the best fit between the ED and the experimental results in the 
region above Tp. 

Figure 5(a) shows the temperature dependence of the x &!■ the model shown in Fig. 4(a). The 
fitting of the data leads to the spin exchange coupling Jep is 183K, while relatively small ones are 
Jg(^ = Jcp = Jgfi =97K, respectively. The ratio ^ed/^ep is 0.53. The (^-factor is 1.71, which is 
rather small compared to (7 ~ 2.0. In principle, the mechanism which reduces the value of the g 
factor estimated from x i^i experiments could exist. One possibility is an effect of small spin-orbit 
couplings. However, we do not discuss this point in this paper. As a reference, the data not only 
for 16-site system but also for 8-site system are shown in Fig. 5(a). The size dependence appears 
below Tp. However, the finite size eff'ect is small above Tp, which justifies the present estimate of 



the spin exchange couphngs to take as the value in the thermodynamic hmit. 

The spin gap for 8-site, 16-site and 24-site systems were calculated. The size dependence of 
the spin gap is shown in the inset in Fig. 5 (a). If we assume that the fitting function for two- 
dimensional spin gap system has a form as As(A^) = As(oo) + jj, the spin gap extrapolated to 
the thermodynamic limit is estimated to be about lOlK, namely Ag/Jcp = 0.550. The reason for 
taking this fitting function is as follows. The energy dispersion of the triplet excitation near the 
wavenumber with the lowest energy excitation is approximated as eT(^) = ^s + c^fe^ + Cyky for 
sufficiently small kx and ky, where Ag represents the bulk spin gap. Using the periodic boundary 
conditions, the wavenumber kx{ky) has a L~^{L~^) correction, where Lx{Ly) is the number of sites 
in the x(y)-direction. Because L^ and Ly are taken as N^^'"^ in two dimensions, the above fitting 
function is obtained. If the energy dispersion can be approximated by the quadratic form in a wide 
range of the wavenumber space, the fitting function reduced from the quadratic dispersion may be 
available even in small size systems. This model may belong to this case. 



Fig. 5. Temperature dependence of uniform magnetic susceptibility for five models. Figures (a), (b), (c), (d) and 
(e) correspond to x for the models in Figs. 4 (a), (b), (c), (d) and (e), respectively. Bold-gray lines represent the 
experimental result. Insets show the extrapolation of the spin gap as a function of A'^ in the ED results. 



For the case shown in Fig. 5(b), we take the relevant spin exchange couplings as J^^^ > Jep, 
for example, Jep = O.OJg^^. The reason is discussed in the next subsection. Figure 5(b) shows 
the temperature dependence of the x- In ^ similar way, the relevant spin exchange couplings are 
estimated to be J^^ ^ 185K and Jep ~ 166K, while the relatively small ones are Jep = Jed = 85K. 
The ratio Jep/ Jed i^ 0.46. The g-factor is 1.72. In this case, one should be careful in estimating the 
spin gap from the data only for small size systems since the system seems to be near the critical 
point where the spin gap vanishes. Beyond the critical point when the spin gap disappears, the 
energy dispersion of the triplet states near the lowest excitation is expected to have a form as 



e'p(fc) = J Cxkx + Cyky, namely, a linear dispersion. As the system with the spin gap is close to the 
critical point, the wavenumber space where the form of the energy dispersion is approximated to 
be quadratic becomes narrower and narrower while the linearly dispersive region becomes wider. 
Then the leading term of the finite-size correction may crossover from N^"^'"^ for small system size 
to N^^ for large system size. Here we take the form As{N) = As(oo) + -j= as the fitting function. 
We have tried extrapolations by using two types of the fitting function. However, we found that the 
leading term of the finite size correction is iV~^' ^ rather than N^^ in this case due to the limitation 
of the system size. Using this fitting function, the spin gap extrapolated to the thermodynamic 
limit is estimated to be about 34K, or Ag/J^^ = 0.186. From the above reason, this estimate of 



the spin gap is expected to give a slight underestimate. In this paper, in order to estimate the bulk 
spin gap, these two fitting functions are employed according as the situation. 

In Fig. 4 (c), the relevant spin exchange couplings are divided into two groups. One is the nearest- 
neighbor spin exchange coupling J and the other is the next-nearest-neighbor one J' . We consider 
the case that J' is larger than J and the other relatively small exchange couplings are neglected. 
The reason is discussed in the next subsection. Then the parameters J and J' are estimated to be 
90K and 237K, respectively. The ratio J / J' is 0.38. The g-factor is 1.68. The relatively large spin 
gap survives in the thermodynamic limit and is estimated to be Ag/J' — 0.564 or 133K when we 
use the fitting function As(A^) = As(cxd) -|- -^. 

The values of the spin gap for the first and the third model are close to the one in the experimental 
results. However, in the second case, the spin gap value is much smaller than the experimental one. 

3. 2 Origin of Spin Gap and Triplet Energy Dispersion 

We consider the mechanism of the spin gap formation. To investigate the dispersion of triplet 
excitations, the perturbation method is useful since it is found by the analysis mentioned in the 
previous subsection that the strength of the relevant spin exchange couplings is relatively large 
compared to the others. In the case shown in Fig. 4 (a), the term in eq.( |2.lD including the relevant 
spin exchange coupling Jep is treated as the unperturbed Hamiltonian, while the other terms are the 
perturbed ones. The strength of the other spin exchange couplings is taken as J . The singlet ground 
state of the isolated edge-shared plaquette is the resonating valence bond (RVB) singlet. Then the 
ground state of the unperturbed Hamiltonian is the product state of the RVB singlets on the edge- 
shared plaquettes. The origin of the spin gap is understood by the edge-shared plaquette singlet. 
The lowest excited state of the plaquette is the extended state of the triplet pair in the RVB state. 
Then the first excited states of the unperturbed Hamiltonian are constructed from the triplet state 
on one of the plaquettes and the singlets on the others. The degeneracy is lifted by the first-order 
perturbation through the transfer of the triplet due to the translational invariance, which yields 
the dispersion of the triplet states. We calculate the energy difference between the ground state 
and the triplet states within the second-order PE. The detailed procedure is explained in Appendix 
A. Figure 6 (a) shows the energy dispersion of the triplet excitation As(fc). The wavenumber 
shown as {kx,ky) represents the vectors kxa2 + kyb2 with 0.2 = cii + 61 and &2 = —(^i + ^i- 
The first magnetic Brillouin zone is expanded by a2 and 62 as shown in Fig. 2. In this case, the 
wavenumber of the lowest excitation is (0,0). Since a strong antiferromagnetic correlation due to 
Jcp is larger than the one due to Jg^^ in this model, the triplet states which have a strong in-phase 
correlation between the nearest-neighbor edge-shared plaquettes become the lowest excited states. 
The amplitude of the spin gap at (0,0) calculated by the second-order PE is 0.575Jg(j, which is 
close to the one estimated by the ED and is also consistent with the experimental result. These 
behaviors qualitatively reproduce the experimental results. 



Fig. 6. Wavenumber dependence of triplet and singlet excited state energies for three models. Figures (a), (b) and 
(c) correspond to the models shown in Figs. 4 (a), (b) and (c), respectively. Filled circles with error bar represent 
the experimental results.cJ^ 



In the case shown m Fig. 4 (b), the terms with the relevant spin exchange couphngs Jep and 
Jp^l are treated as the unperturbed Hamiltonians, while the other terms are the perturbed ones. 
In the Jg(| > Jep case, the product state of the dimer pairs on J^^ bonds is the unperturbed 
ground state. Then the origin of the spin gap is ascribed to the corner-shared dimer singlets. The 
lowest excited states of the edge-shared plaquette are the product states of the triplet pair on one 
Jpj^ bond and the singlet pair on the other J^^ bonds. Because two kinds of triplet states exist, 
the lowest excited states are six-fold degenerate, which is also different from the previous case. We 
calculate the energy gap between the singlet state and those triplet states in a similar way. Figure 6 
(b) shows the energy dispersion of the triplet states. The degeneracy is lifted by the first-order 
perturbation and the bonding and the anti-bonding triplet states appear. The wavenumbers of 
the lowest excitations are (0, vr) and ('7r,0). The spin gap is 0.197 J^^j^, which is consistent with 
the value obtained by the ED. In this case, the square area surrounded by the lines connecting 
the nearest-neighbor points of the lowest excitation becomes a half of the original first magnetic 
Brillouin zone in the wavenumber space, which contradicts the experimental results. 

In the model shown in Fig. 4 (c), the ground state of the Hamiltonian given by the relevant spin 
exchange couplings is not straightforwardly obtained since the system has no isolated terms. To 
understand the mechanism of the spin gap formation, we neglect the small spin exchange couplings 
other than J and J' in this model. We consider three possible ground states and the origins 
of the spin gap for the Hamiltonian with only the relevant spin exchange couplings. To clarify 
three possibilities with the help of the perturbation method, the relevant terms are divided into 
two parts, unperturbed and perturbed terms. In the first case, the unperturbed Hamiltonian is 
taken as the term with Jep while the others are the perturbed ones. Then the ground state of 
the unperturbed Hamiltonian is represented by the edge-shared plaquette singlet. In the second 
case, we take the unperturbed Hamiltonian as the Jg^j term. Then the ground state is described 
by the edge-shared dimer singlet. In the third case, the unperturbed Hamiltonian is the term with 
the next-nearest-neighbor spin exchange couplings. Then the ground state is represented by the 
four-site stripe singlet. Since the first and the second cases have been already investigatedpEP we 
concentrate on the third case. Then the amplitudes of the nearest and the next-nearest neighbor 
couplings are chosen as J and J', respectively. Figure 7 shows the ground state energy per site 



as a function of the strength of the spin exchange couphng J /{J + J') calculated by the PE. The 
detailed calculation by the PE is written in Appendix B.The ground state described by the stripe 



Fig. 7. Ground state energy per site eg/(J + J') as a function of J /{J + J'). Bold, broken and dash-dotted lines 
represent the ground state energy by the second-order perturbation from the plaquette singlet, dimer singlet and 
stripe singlet, respectively. 



singlet is favored when J /{J + J') is smaller than 0.483, or the ratio J/ J' is smaller than 0.932. 
The system with J = 90K and J' = 237K estimated by the ED belongs to the stripe singlet phase 
as shown in Fig. 7. Therefore, we may conclude that the stripe singlet is a possible mechanism of 
the spin gap formation in this model. We calculate the triplet dispersion as shown in Fig. 6(c). 
The ground state is the product state of the stripe singlet for the unperturbed Hamiltonian. The 
magnetic unit cell in the real space contains two stripes crossing each other. Then the lowest 
triplet states in the unit cell are six-fold degenerate. Similarly to the analysis for the model given 
in Fig. 4(b), the bonding and anti-bonding triplet states appear. The spin gap is 0.539 J', which 
is close to the value obtained by the ED and is also close to the experimental result. However, 
the wavenumbers of the lowest excitations are vrcii and vrbi respectively, which is inconsistent with 
the experimental results. The magnetic periodicity is also different from the experimental one, 
since the magnetic first Brillouin zone is the same as that of the unit lattice. In this model, small 
exchange couplings in the corner-sharing have been neglected. They may change the quantitative 
feature such as band-width of the triplet dispersion. However, an effect due to them is very small 
since both of the occupied d orbitals on V atoms do not extend to the oxygen. Qualitative aspects 
such as the origin of the spin gap and the magnetic periodicity are not affected by these small 
exchange couplings. 

3.3 Scattering Intensity 

We investigate the scattering intensity of the neutron inelastic scattering given by the experi- 
ments. The scattering intensity is proportional to the Fourier component of the spin-spin correlation 
written as 

I(q,a;) ^ldte''^\Siq,tyS{-q,0)), (3.2) 

where 

5(q,t)=iV-^^5,(t)e^i9-n^ (3.3) 

i 

and (• • •) represents the thermal average. The wavenumber q are expanded by the vectors di and 
bi, namely we take Qx axis in the di direction and qy in the bi direction. When the temperature 
is sufficiently lower than the spin gap, eq.(|3.2D is approximated by the transition probability from 
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the ground state |g.s.) to the triplet states |n), which is described as 

I{q,u:) oc ^ |(n|5(q)|g.s.)|25(w - S„ + Eg), (3.4) 

n 

where the summation is restricted to the triplet states and En and Eg are the energy of the triplet 
states and the ground state, respectively. Here, when we calculate eq.( |3.4| ), we use the wavefunctions 
of \n) and |g.s.) obtained by the PE, which is written in Appendix A in detail. The integrated 
intensity of the neutron inelastic scattering is obtained theoretically by integrating eq.(|3.4|) over uj, 
described as 

I[q)^ f dwl{q,u:) ^Y.MS{q)\g.s.)\\ (3.5) 

Figures 8(a), (b) and (c) show the wavenumber dependence of the integrated intensity calculated 
by the PE. The experimental results are also shown in these figures. Especially, Figs. 8(a) and (c) 
indicate important points. In the case shown in Fig. 4(a), since the nearest-neighbor antiferromag- 
netic correlation due to Jep is dominant, the scattering intensity has a minimum at {'K/a)ex and 
becomes larger as the wavenumber approaches {7r/a)ex + {■K/a)ey. While in the cases shown in 
Figs. 4(b) and (c), since the next-nearest-neighbor antiferromagnetic correlation is dominant, the 
scattering intensity has a maximum at (7r/a)e^ and becomes smaller as the wavenumber approaches 
[t: /a)e,j. + {■K/a)ey. From this result, experimental results indicate that the next-nearest-neighbor 
antiferromagnetic correlation is strong.ll2P This supports the models in Figs. 4(b) and (c). How- 
ever, the wavenumber dependence for the model in Figs. 4(c) is qualitatively different from the 
experimental results as shown in Figs. 8(a) and (b). Therefore, the experimental results are rather 
consistent with the integrated intensity for the model shown in Fig. 4(b) than those for the models 
shown in Figs. 4(a) or (c). 



Fig. 8. Wavenumber dependence of the integrated intensity for three models shown in Fig. 4(a), (b) and (c). The 
wavenumber q in Figures (a), (b) and (c) are represented by q = 27r(ai -\-kbi) with 0.5 < k < 1.5, q — 2n{hdi +bi) 
with 1 < h < 2 and q = 27r[(l + 5)di + (1 — ^)''i] with < 5 < 1, respectively. Filled circles represent the 
experimental results given by Kodama et alx3 The right ordinates represent the scale of the integrated intensity of 
the neutron inelastic scattering experiments. Here, the value 200 counts/5500kmon in the experiments corresponds 
to 0.04 obtained by the PE. This gives the best fit of the data for the model shown in Fig. 4(b) with the experimental 
results. 



§4. Discussion 

To understand the mechanism of spin gap in CaV409, we have studied the orbital-order depen- 
dence of the physical quantities such as the uniform magnetic susceptibility, the energy dispersion 
of the triplet states and the integrated scattering intensity. Before discussing these results, we show 
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the temperature dependence of x for the system with Jcp = 171K, J^^^ = 15K and Jep = Jg^j = 68K. 
These parameter values are the ones suggested by an analysis of the neutron inelastic scattering 
experiments.tJ^ The uniform magnetic susceptibilities for 8-site and 16-site systems are calculated 



Fig. 9. Temperature dependence of the uniform magnetic susceptibility for the model with Jcp = 171K, J„j — 15K 
and Jep — Jp^ — 68K. Broken and solid curves represent x for 8-site and 16-site systems calculated by the ED, 
respectively, while the bold-gray line is the experimental result. 



by the ED. The theoretical and the experimental results are shown in Fig. 9. The size dependence 
is almost negligible and the numerical result of the x for 16-site system may be regarded as the 
one in the thermodynamic limit. Since it is difficult to estimate the amplitude of the spin exchange 
couplings and the (7- factor from the high temperature expansion, the g-factor is taken as 1.7 here, 
which is nearly the same as the value obtained by the analysis mentioned in the previous section. 
Figure 9 indicates that a qualitative difference between the experimental and theoretical results 
appears around the peak temperature. We may conclude that the parameters suggested by the 
analysis of the neutron inelastic scatteringl^P are not consistent with the experimental result of the 
uniform magnetic susceptibility. 

In Section 3, we have mainly studied the cases in Figs.4(a)-(c). Here, we briefly discuss the cases 
in Figs. 4(d) and (e) for completeness. Figure 5(d) and (e) show the temperature dependences of x- 
For the model in Fig. 4(d), strength of the relevant spin exchange coupling Jj^ and the other small 
ones J are estimated to be 159K and lllK, respectively. The ratio J/Jji is 0.70. The /7-factor is 
1.71. The amplitude of the spin gap is extrapolated to be about 29K or As/Jr ~ 0.179. For the 
model in Fig. 4(e), Jj^ and J are estimated to be 207K and 89K, respectively. The ratio J/Jji is 
0.43. The g-factor is 1.70. The amplitude of the spin gap is about 3.6K or As/Jr ~ 0.017. In 
both cases, the model with given parameters Jji and J seems to be near the critical point because 
the leading term of the finite size correction in the fitting function is N^^'"^ rather than N~'^ and 
the extrapolated values of the spin gap are very small. From the perturbational results, the origin 
of the spin gap for these models may also be described by the stripe singlet as well as the case in 
Fig.4(c). 

Here, we discuss more quantitative aspects of the magnetic properties. We first discuss the 
temperature dependence of x- The experimental results for the x have not been well understood 
by the model with only the nearest-neighbor interaction, i.e. Jcp = J^fj = 0. In experiments, the 
spin gap As, the peak temperature Tp and the temperature T* which gives a half amplitude of 
X at Tp are 107K, llOK and 595K, respectively. The ratios Ag/Tp and T* /Tp are 0.97 and 5.41, 
respectively. The peak temperature is relatively small compared to the spin gap and x decays 



12 



rapidly below Tp. The numerical results for the model with the nearest-neighbor interactions have 
several different features from the experimental ones.BEP For example, the ratios Ag/Tp and T* /Tp 
for the isolated plaquette model(Jg(^ = 0) are 1.28 and 3.83, respectively. In the plaquette singlet 
regions, these ratios become small as Jg^j becomes large. While the ratios Ag/Tp and T* /Tp for the 
isolated dimer model (Jep = 0) are 1.60 and 3.48, respectively. In the dimer singlet regions, these 
ratios also become small as Jep becomes large. Then the ratios estimated from the experiments are 
no longer reproduced from this model. We consider several possible reasons of this discrepancy from 
the experimental results below Tp. One is that the frustration due to the next-nearest-neighbor 
couplings reduces the band width of the triplet dispersion while keeps the strength of the spin 
gap, as has already been pointed out.EP The other is that the low- lying excitations from the singlet 
ground state are not only the triplet states but also the singlet states. For example, in the case 
shown in Fig. 4(b), the singlet dispersion exists in the low energy region, as shown in Fig. 6(b). 
These singlet excited states are constructed from the product states of the lowest excited singlet 
state on one of the edge-shared plaquettes and the corner-shared dimer singlets on the others for 
the unperturbed Hamiltonian. Within the PE, the energy gap between the lowest singlet excitation 
and the ground state is comparable to the spin gap. In addition, since the singlet dispersion is 
nearly flat, the energies of the singlet excitations are lower than those of the triplet excitations in a 
wide region in the Brillouin zone. Then at the temperature around the singlet-singlet gap, x decays 
rapidly since the strength of the denominator in eq. ( |3.lD increases relatively due to the increase of 
the weight of the excited singlet states. In order to discuss the possibility of the existence of the 
low-lying excitation of the singlet states, the detailed analysis of the experimental results for the 
specific heat is needed because the temperature dependence of x below the spin gap temperature 
or the neutron inelastic scattering do not contain the information about the singlet excitations 
directly. 

Next, we discuss the treatment of the spin exchange couplings. The spin exchange couplings 
have been divided into two groups. One is strong and the others are relatively small. The same 
value has been assigned for the values of the other relatively small couplings. This weaker bonds 
have minor contributions to physical properties and do not change the essential feature obtained 
in Sec. 3. However, within these weaker bonds, it is possible that the amplitude of the nearest- 
neighbor couplings becomes twice as large as that of the next-nearest-neighbor couplings. Because 
the number of the path for the superexchange mechanism is two in the nearest-neighbor case while 
it is one in the next-nearest-neighbor one. If this estimate is adopted, for example in the case 
shown in Fig. 4(a), the band width of the triplet dispersion becomes narrow and the wavenumber 
giving the lowest excitation shifts from (0,0) to (tt, vr) within the second-order PE. Since, in the 
terminology of the perturbation method, the transfer energy of the triplet between the nearest- 
neighbor edge-shared plaquettes is determined by Jg^j — 2 Jep within the first-order PE, it vanishes 
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at Jgj^ = 2Jcp. The higher-order perturbation may change the wavenumber of the lowest excitation 
according as the sign of the transfer energy of the triplet or the amplitudes of the nearest and the 
next-nearest neighbor couplings. 

Through this study, we have concentrated on the superexchange mechanism for the dxz and the 
dyz orbitals. In the real materials, however, there are several possible subtleties which further 
invalidate the assumption that the amplitude of the relatively small nearest-neighbor couplings is 
twice as large as that of the relatively small next-nearest ones. One is the contribution from the 
superexchange mechanism in the d^y orbital since the wavefunction of the occupied orbital may also 
have a dxy component in the real material. Due to this contribution, the next-nearest-neighbor 
couplings become antiferromagnetic and relevant. Another is the contribution from the direct 
exchange mechanism for both the dxy orbital and the dxz and dyz orbitals. Since the bonding 
of the t2g orbitals with the pz orbital on oxygen is not a bonding but tt bonding, the effect of 
the direct exchange mechanism may be relatively strong compared to that of the superexchange 
mechanism. The detailed character of the direct exchange couplings may depend on the symmetry 
of the wavefunction and the distance between the V atoms. The contributions from tilting of 
the pyramid or the Jahn- Teller distortion may also influence the exchange coupling. Detailed 
quantitative analysis for the spin exchange couplings remain for further studies. 

Finally, we discuss the importance of the dxz and dyz orbitals in the magnetic properties for 
CaV409. In the real material, the wavefunction lying near the Fermi energy level may be the 
linear combination of the d orbitals. In this context, Marini and Khomskiit^^ also pointed out 
the importance of the orbital order. They discussed the effect of the crystal field and estimated 
the coefficient of the linear combination of the dxz and dyz orbitals in the wavefunction. The 
given stable orbital is tilted and has a chirality. They concluded that the relevant spin exchange 
coupling is Jg^j and the origin of the spin gap is an edge-shared dimer singlet. They proposed that 
the next strongest spin exchange couplings are the next-nearest-neighbor ones and Jep is weakly 
ferromagnetic. In this case, however, the nearest-neighbor spin correlation is enhanced and the out- 
of-phase correlation between the nearest-neighbor edge-shared dimer singlets should be strong due 
to the next-nearest-neighbor antiferromagnetic spin exchange couplings and weakly ferromagnetic 
Jep. These contradict the reported experimental results.c3 

Neutron measurement suggests substantial antiferromagnetic correlations for the next-nearest- 
neighbor pair of V. This is naturally explained by the occupation of dxy orbital. However, if only 
dxy orbitals are occupied, it is hard to justify the model by Kodama et alS3 where a much larger 
Jep than Jg(j is assumed. Aside from details it appears indeed necessary to take different Jep 
and Jg(j to reproduce the nuetron data. Our model of orbital order can explain why the two 
next-nearest-neighbor spin exchange couplings differ. However, simple model with only dxz and 
dyz occupations appears to be insufficient to explain both of the susceptibility and neutron results 
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simultaneously as we discussed in this paper. All the above results imply, within this framework, 
that the wavefunction is in fact represented by a linear combination of the dxy orbital and the 
other t2g orbitals with the orbital order. We propose that the ground state is described as the 
wavefunction with large weight of dxy orbital superposed with the ordered dxz and dyz orbitals 
given in Fig. 4(a). In this case, the basic origin of the spin gap may be a corner-shared plaquette 
singlet. Then the next-nearest-neighbor antiferromagnetic spin correlation becomes relatively large 
since the large weight of dxy orbital enhances Jcp, which may explain the experimental results of 
the wavenumber dependence of the integrated intensity of the neutron scattering. Additionally, a 
strong in-phase correlation between the nearest-neighbor corner-shared plaquette singlets increases 
because of the existence of Jg^j due to the presence of the orbital order. Then the wavenumber of 
the lowest triplets may become (0,0), which is consistent with the recent experimental results.Ej^ 
If this situation is true, small plaquette which is constructed by four oxygens without a V atom at 
the center may rotate or some other specific type of the lattice distortion may occur in order to 
favor the given orbital order. To obtain the information on the weight of these orbitals, a detailed 
analysis of the lattice structure is needed, especially at very low temperatures. Quantitative and 
detailed theoretical analyses of this proposal remain for further studies. 

The orbital order effect may also be observed in the temperature dependence of x for CaV205.o) 
This material also shows the spin gap behavior. The spin gap and the peak temperature are about 
464K and 400K, respectively. One possible mechanism of the spin gap formation may be due to 
the ladder structure. Since it is assumed that d electron is localized at the dxy orbital, the spin 
exchange couplings J in the leg and J± in the rung have almost the same values. Numerical 
analysis of x foi' the AFH model on the ladder lattice has led to the result that the ratio As/Tp is 
almost 0.5 at J = Jj_.li3 However, the ratio obtained in the experiments is 1.16, which is different 
from the theoretical prediction. If the localized orbital on V atom is represented by the linear 
combination of the dxy and one of the other t2g orbitals which extends to the oxygen on the rung 
due to the orbital order, J± becomes larger than J. Since the ratio Ag/Tp is 1.60 in the dimer 
limit (J = 0), it may become 1.16 as J increases from zero. Quantum or thermal fluctuation of 
the orbital degree of freedom also changes a quantitative behavior of x at finite temperature. We 
propose that the orbital order of dxz and dyz orbitals with partially fillings of t2g orbitals is also a 
promising explanation for rather puzzling experimental results in CaV205. The quantitative study 
for this problem remains for further study. 

§5. Summary 

We have investigated orbital order effects in magnetic properties of CaV409. Several possible 
models with the orbital order have been considered. The amplitudes of the spin exchange couplings 
for each model are determined in order to reproduce the temperature dependence of x ^^ experi- 
ments. The origins of the spin gap are not necessarily the same as the originally proposed plaquette 
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singlet but are ascribed to the generalized four-site singlet. By using the estimated values of the 
spin exchange couplings, the dispersion of the triplet excitations and the spin-spin correlation cor- 
responding to the integrated scattering intensity have been calculated within the PE. The strong 
in-phase correlation between the nearest-neighbor pair of four-site singlets and the large antifer- 
romagnetic spin correlation between the next nearest-neighbor spins are minimal requirements for 
the explanations of the experimental results in the neutron scattering. Although the wavenumber 
dependences of the triplet dispersion and the spin-spin correlation for our models do not show com- 
pletely satisfactory agreement with those obtained by the presently available experimental data, 
they are much improved from the single-orbital cases. We further propose the mechanism required 
to explain each feature of the experimental results. We suggest that the order of the dxz and dyz 
orbitals hybridized with uniform and partial occupation of dxy orbital may explain the experimental 
results. We have also discussed that effects of the orbital order may be observed in other vanadium 
oxide compounds such as CaV205. These strongly suggest that the role of the orbital order is 
important in understanding the magnetic properties of the vanadium oxide compounds with the 
spin gap. 
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Appendix A: Triplet Dispersion and Wavefunction by Perturbation Expansion 

We introduce the procedure for calculating the energy dispersion of the lowest triplet band and 
the wavefunctions of the singlet ground state and the triplet band. 

We first consider the model shown in Fig. 4(a). The unperturbed Hamiltonian has the terms with 
Jep on the edge-shared plaquette, while the other terms are treated as the perturbed ones. The 
ground state of the edge-shared plaquette with only Jep is the RVB state \(j)g). The energy of this 
state is —2 Jep. Then the ground state of the unperturbed Hamiltonian is the product state of \(j)g) 
written as 

l^g)(a)=ni'^g(^^))- (A-1) 

i 

The lowest triplet states of the edge-shared plaquette \ip) with Sz = —1, 0, 1 are the hopping state 
of the triplet pair in the RVB state. The energy of these states is —Jep- Then the wavefunctions 
of the lowest triplet band are constructed by \ip) on one of the edge-shared plaquette and \4>g) on 
the others. The degeneracy is lifted by the first-order perturbation. The wavefunctions of these 
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triplet states are written as 



|^(fc))(,) = -^ E \^^^i)) n \<t>g{r,))e'^ ■ ^% (A.2) 



where Np is the number of the edge-shared plaquette. Within the second-order PE, the energy 
difference between the triplet states and the singlet ground state is represented as 

J2 j2 

(cos 2kx + cos 2ky) — -^^ — cos kx cos ky. 



108 Jep^ '' 9Jep 

Next, we consider the model shown in Fig. 4(b). The unperturbed Hamiltonian has the terms 
with Jep and J^^^ on the edge-shared plaquette while the other terms are the perturbed ones. Since 
we take J^.^^ > Jep, for example, Jep = 0.9J(,(^, the ground state of the edge-shared plaquette is 
not the \(j)g) but the product state |^g) of the dimer pair on each J^^^ bond. The energy of this 
state is —1.5J^^. Then the ground state of the unperturbed Hamiltonian is the product state of 
|,^g) written as 

i^g)w=ni^g(^^))- (A-3) 

i 

The wavefunctions of the lowest triplet states on the edge-shared plaquette are constructed by two 
states \ri) and |0 with Sz = —1,0, 1 which are the product states of the triplet pair on one of the 
Jpj bond and the dimer pair on the others. Then the wavefunction of the lowest triplet bands are 
constructed from one of the triplets \rj) or \0 on one of the edge-shared plaquette and |^g) on the 
others. The first-order perturbation lifts the degeneracy. The bonding and anti-bonding states of 
the triplets appear, which are written as 

|^I/±(fc))(,) = -^ Y.(^^^\vir^)) + 7^^"^ \C{n))) n ieg(r-,))e^^ ' ^% (A-4) 

V^^P k V10T4V5 V10T4V5 j^i 



in the parameter regions < k^ < ky < it. Here Jep and J^j are both taken as J. The energy 
difference between these triplet states and the ground state is obtained as 

3J2 J2 3J2 5J2 



As(fc)± = Jed + 



+ ^ + 



-1 =F Vo) cos kx + (—1 ± v5) cos k 



8Jcd 16(Jcd - «^ep) 8(2Jed " <^cp) 8(2Jed + Jep) 
'J ,P ' 

4 + 8Jed, 

In this model, the lowest singlet band is also calculated. In a similar way, the energy difference 
between the lowest singlet band and the singlet ground state is obtained as 

9J2 3J2 J2 j2 

A's(fc) = 2( Jed - ^^cp) + ^^ - ^^ - ^.j , J . + ^l-(cos kx + cos ky). (A-5) 

»Jed o<^ep ^Wcd + "^epj oJep 
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Next we introduce the wavefunctions for the models shown in Fig. 4(c). In the case shown in 
Fig. 4(c), the unperturbed Hamiltonians are taken as the part with the relevant next-nearest- 
neighbor spin exchange couplings, while the other part with the relevant nearest-neighbor ones 
are the perturbed ones because we suggest that another possible origin of the spin gap is described 
as the stripe singlet in this case. Since the magnetic unit cell has two stripes, the ground state 
of the unperturbed Hamiltonians |<l>g)(c) is represented by the product state of the stripe singlet 
Wg-iTi)) 1,2 on each stripe written as 

\^g)(c)=\{\4>'g{r^))l\4>'g{r.^)2. (A-6) 

i 

where the indices 1 and 2 represent the one stripe and the other stripe in the magnetic unit cell, 
respectively. The lowest triplet states of the magnetic unit cell for the unperturbed Hamiltonian are 
six-fold degenerated i.e. |V'')il'^g)2 and |i?i>g)i|V'')2 with Sz = —1,0, 1. Then the bonding and anti- 
bonding triplet bands appear by the PE. Within the first-order perturbation, the wavefunctions 
are obtained as 

|vI/±(fc))(,) = -^5:(|^'(r,))i|<A'g(n))2±|0^(r,))i|V''(n))2)ni0k(^^-))il0g(^.))2,ei^-^% 

(A-7) 
where A^s is the number of the stripe. We calculate the energy difference between these triplet 
bands and the ground state within the second-order PE written as 

As(fc)± = J'[Ai ± A2{- cos kx + cos ky)], (A-8) 

where the parameters Ai and A2 are 0.601 and 0.031, respectively at J / J' = 0.38. Here the 
component kx and ky are represented by the vectors di and &i, respectively. 

Appendix B: Ground State Energy for Model in Fig. 4(c) 

We calculate the ground state energy per site by the second-order PE to investigate the possible 
mechanism of the spin gap formation for the model shown in Fig. 4(c). The first case is that the 
unperturbed Hamiltonian is the term with Jep and the other terms are the perturbed Hamiltonian. 
In this case, the origin of the spin gap for the unperturbed Hamiltonian is the edge-shared plaquette 
singlet. The strength of the other spin exchange couplings is taken as J. The ground state energy 

per site is obtained as 

1 1 11 J2 , ^ 

e„ = — J +—J . B-1) 

S 2 °P 16 576 Jep ^ ' 

The second case is that the unperturbed Hamiltonian is the term with J^^j while the others are the 

perturbed ones. The strength of the other spin exchange couplings is taken as J. In this case, the 

origin of the spin gap is the dimer singlet. The ground state energy per site is obtained as 

3 9/2 

e„ = --J , - — — . (B-2) 

^ 8 ^^ 128 Jed 
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The third case is that the term with the next-nearest-neighbor spin exchange couphng J' is the 
unperturbed Hamiltonian while the others are the perturbed ones. Similarly, the strength of the 
other spin exchange coupling is taken as J . In this case, the origin of the spin gap is the stripe 
singlet. The ground state energy per site is obtained as 



^ 16 J' V8 4 

where the parameters A and B are written as 



768 



A = —, (B-4) 

3(1 + ^/3)' ^ ' 



(4 - ^/2 - 2^/3 -F \/6)2 (2-3^/2 + ^/6)2 8(^/3-1)2(2-^/2) 
l-^/2 + ^/^ l + \/3 2-^/2 + 2^/3 



(-4-^2 + 2^3 + ^6)2 (2 + 3^/2 + ^6)2 8(^3-1)2(2 + ^2) 
l + \/2 + V3 l + \/3 2 + ^/2 + 2^/3 

respectively. When the parameters Jep and Jg^j are taken as J, the J /{J + J') dependence of the 
ground state energy per site is obtained, as shown in Fig. 7. 
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